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Some new developments obtained in the last few years concerning the propagation of high energy cosmic rays are 
discussed. In particular, it is shown how the inclusion of drift effects in the transport diffusion equations leads 
naturally to an explanation for the knee, for the second knee and for the observed behavior of the composition and 
anisotropics between the knee and the ankle. It is shown that the trend towards a heavier composition above the knee 
has significant impact on the predicted neutrino fluxes above 10 14 eV. The effects of magnetic lensing on the cosmic 
rays with energies above the ankle are also discussed, analyzing the main features of the different regimes that appear 
between the diffusive behavior that takes place at lower energies and the regime of small deflections present at the 
highest ones. 



1 Introduction 

Since their discovery in 1912, cosmic rays (CRs) were 
of great help for particle physics, providing a source 
of high energy particles for free, which only required 
the construction of detectors in order to observe dif- 
ferent kinds of interesting phenomena. In this way, 
positrons, muons, pions, kaons and hyperons were 
discovered in the period 1930-1950. However, when 
after the '50s man made accelerators reached ener- 
gies beyond the GeV, particle physics moved back 
to the labs and cosmic ray research became focused 
on the study of the CRs themselves (rather than on 
the products of their interactions), trying to under- 
stand their origin, the mechanisms responsible for 
their acceleration and the way they propagate from 
the sources up to us. 

There are only a few observable quantities asso- 
ciated to the CR fluxes. These are the energy spec- 
trum, the CR composition and the anisotropies in 
arrival directions, and it is through their study that 
the CR mysteries have to be unraveled. For instance, 
the fact that the spectrum is essentially a power law 
and not a thermal one is what led Fermi to suggest 
that the CR acceleration was a stochastic process. 
Also, looking at low energy cosmic rays, the study 
of the abundances of spallation products (like Li, Be 
and B, which are produced mainly by spallation of 
C, N and O) gives information about the amount of 
matter traversed by the CRs, while the abundances 
of unstable isotopes (e.g. 10 Be) gives information on 
the time spent by CRs in the Galaxy. On the other 
hand, the anisotropies observed on the low energy 
CRs arriving from the East and from the West gave 



indications that they were caused by the deflections 
produced by the geomagnetic field on the positively 
charged CRs. 

Many puzzles are also associated with the ob- 
served properties of the CRs with very high ener- 
gies, those above 10 15 eV and up to the highest 
ones exceeding 10 20 eV. In particular, we would like 
to know what causes the spectral changes observed, 
what is the origin of the observed anisotropies and 
the changes in composition as a function of energy 
and how CR sources would look like at ultrahigh en- 
ergies. As will be discussed below, a crucial issue in 
this respect is to understand in detail the propaga- 
tion of CRs through the magnetic fields present in 
the Galaxy, since this is essential to determine their 
properties when we finally observe them. 

2 The cosmic ray spectrum 

The differential flux of CRs changes by more than 
30 orders of magnitude in the energy range from 
10 9 eV to 10 20 eV, following essentially a power 
law E~ a which shows only some small but notice- 
able breaks in the power index a. One has indeed 
a ~ 2.6 -7- 2.7 at low energies, with a steepening to 
a ~ 3 at the so-called knee, which takes place at 
Ek — 3 x 10 15 eV. A second steepening to a ~ 3.3, 
referred to as the second knee, has been reported at 
E s k — 4 x 10 17 eV, while a hardening to a ~ 2.7 
takes place at E a ~ 5 x 10 18 eV, being known as the 
ankle. These general features are apparent in the 
data 1 shown in Fig. 1, which displays the quantity 
E 3 x d$/dE. The observed spectrum is just a reflec- 
tion of the original source spectra which have been 
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Figure 1. CR spectrum as measured by different experiments 
and predictions of the diffusion/drift scenario for different 
galactic models. The dotted line is the assumed extragalactic 
flux. 



reshaped by the energy dependence of the confine- 
ment time of CRs in the Galaxy, since for increasing 
energies CRs escape more readily and hence their lo- 
cal density becomes less enhanced. To understand 
this last effect we will then discuss in some detail the 
transport of CRs through the magnetic fields perme- 
ating the Galaxy. Beyond the ankle the CRs arc no 
longer confined inside the Galaxy (and most likely 
there are no galactic sources reaching such high en- 
ergies) so that a transition to an extragalactic com- 
ponent should be taking place there. 

2.1 Cosmic ray diffusion and drift: 

Our Galaxy is believed to have a regular and a ran- 
dom magnetic field, both with local strength of a few 
yuG. The regular component is aligned with the spi- 
ral arms, reversing its directions between consecutive 
arms and having a typical scale height of ~ 1 kpc. 
Also a more extended halo component is believed 
to exist, with typical scale height of a few kilopar- 
secs but its symmetry properties are less well es- 
tablished. On the other hand, the random compo- 
nent has a spectrum not very different from a Kol- 
mogorov one (as would result if its origin is associ- 
ated with the turbulence in the interstellar medium 
(ISM)), i.e. with a magnetic energy density satisfy- 
ing dE r /dk oc fc~ 5 / 3 in Fourier space, with a maxi- 
mum scale of turbulence of order L max ~ 100 pc. 



When particles of charge Ze propagate across 
a regular field B , they describe helical trajectories 
characterized by a pitch angle 9, so that the compo- 
nent of the velocity parallel to B is v\\ = ccos6*, and 
a Larmor radius given by 



tl = 



pc ^ E/Z ^G_ 
ZeB Q ~ 10 15 cV~Bo~" 



pc 



(1) 



(the radius of the helical trajectory is r^sin^). In 
the presence of the random component B r , the CRs 
will scatter off the magnetic field irregularities with 
associated scales of order tl, changing their pitch 
angle but not their velocity. This will lead to a ran- 
dom walk and a diffusion along the magnetic field 
direction characterized by a diffusion coefficient 



D« = 



(Asjj) 
2At 



.-A,, 



(2) 



where A|| is the pitch angle scattering length, which 
depends on how much power there is in the magnetic 
field modes at scales ~ r^, i.e. 



An oc 



dEr/d Ink 



(3) 



k=2-rr/r L 



Hence, for a Kolmogorov spectrum one has Du oc 
E 1 / 3, , and in general if the spectrum of the random 
magnetic field energy satisfies dE r /dk oc fc m ~ 2 , one 
has £>|| oc E m . 

The diffusion orthogonal to the regular magnetic 
field direction is typically much slower (unless the 
turbulence level is very high, in which case parallel 
and perpendicular motions become similar), and it 
is due to both pitch angle scattering and to the wan- 
dering of the magnetic field lines themselves, which 
drag with them the diffusing particles in the direc- 
tion perpendicular to Bo- Its evaluation has then to 
be done numerically, and recent results 2 show that 
for fixed levels of turbulence (i.e. for given values of 
\J (B 2 ) /Bq) the associated diffusion coefficient D± 
has a similar energy dependence as D\\ as long as 
tl ^ E max . 

The third ingredient is the antisymmetric (or 
Hall) diffusion 3 , which is associated to the drift of 
the CRs moving across the regular magnetic field. 
What is relevant here is the macroscopic drift asso- 
ciated to the gradient of the CR density, which leads 
to a current 



J A = D A ^- x VN. 

-t>0 



(4) 
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Figure 2. Illustration of the physical origin of the macroscopic 
drift. 



Notice that this macroscopic current is orthogonal to 
both Bo and VAT, and its relation to the microscopic 
drift associated to gradients and curvature in Bo is 
subtle 4 . In particular, the macroscopic drift is also 
present for a constant Bo, and its origin is illustrated 
in Fig. 2. The antisymmetric diffusion coefficient is 
just given by 



D A ~ -j- cx E. 



(5) 



The CR density distribution in the Galaxy is then 
determined by the equation 



V • Jd = Q, with J Di 



(6) 



where Q describes the distribution of sources and the 
diffusion tensor is (adopting here the z axis along the 
direction of Bo) 

D ± D A 



Di 



-D A 










Du 



(7) 



To solve equation (6) it is convenient to assume for 
simplicity that the Galaxy has cylindrical symmetry 
and that the magnetic field is azimuthal, what is not 
far from being true, in which case Du plays no role 
and one has 

J D ~ —DxVN + x VN, (8) 

so that the turbulent diffusive component is orthog- 
onal to the isodensity contours while the drifts are 
parallel to them. 

Hence, adopting a realistic configuration of mag- 
netic fields, obtaining from them the diffusion tensor 
everywhere in the Galaxy (using in particular fits to 




Figure 3. CR isodensity contours for E/Z = 10 14 , 3 X 10 15 
and 5 X 10 17 eV (from top to bottom). The arrows indicate 
the directions of the diffusion currents while the shading the 
level of anisotropy. 
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the numerical results obtained for D± for different 
turbulence levels 2 ), assuming a distribution of CR 
sources and then numerically solving the diffusion 
equations, one can obtain 5 ' 6 the CR density distri- 
bution in the Galaxy and the resulting diffusion cur- 
rents. This is exemplified in Fig. 3, which shows the 
isodcnsity contours (every two solid contours corre- 
spond to a change in density by an order of magni- 
tude) for a source of constant strength in the inner 
3 kpc of the galactic plane. The arrows indicate the 
direction of the diffusion currents, with the compo- 
nent parallel to the contours arising from the drifts. 
Notice that since the motion of a charged particle in a 
magnetic field is determined just by its rigidity, the 
relevant variable is the ratio E/Z. The three pan- 
els shown correspond to E/Z = fO 14 , 3 x fO 15 and 
5 x 10 17 eV, and it is clearly seen that for energies be- 
low ZEk the dominant effect arises from the perpen- 
dicular diffusion, while for larger energies the drifts 
are responsible for the dominant escape mechanism. 
The transition between these two regimes is natu- 
rally understood from the different energy depen- 
dence of the two diffusion coefficients (D±_ oc E 1 ^ 3 
while Da oc E). Actually, both coefficients become 
comparable, D± ~ Da, at an energy of the order 
of ZEk, and this crossover naturally generates the 
observed break in the spectrum. 

2.2 The knee: 

If the source spectrum is taken as a power law, 
dQ/dE oc E~ a ° , it is seen from the diffusion equa- 
tion that the observed spectrum will be affected 
by the energy dependence of the diffusion coeffi- 
cients. This just reflects the fact that the CR con- 
finement time in the Galaxy is r e oc D^ 1 . Hence, 
below the knee, where transverse turbulent diffu- 
sion dominates, the observed spectral index will be 
a ~ a s + 1/3, while in the drift dominated regime 
one has instead a ~ a s + 1. Since this transition 
takes place at different energies for CRs with differ- 
ent charges, one expects that first the proton compo- 
nent will steepen its spectrum at an energy that will 
be identified with Ek, and then the heavier compo- 
nents will suffer a similar change in their spectrum 
but at energies ZEk- This is illustrated in Fig. 4, 
which shows the local density of the different CR 
components obtained 5 in the scenario described in 
Fig. 3. The CR composition is normalized with ex- 
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Figure 4. Predicted individual and total spectra in the diffu- 
sion/drift scenario. 

isting satellite measurements of lower energies, typi- 
cally below 10 14 eV, and the individual spectra of H, 
He, O and Fe are displayed, together with the total 
flux. 

ft is seen that although each individual compo- 
nent changes its spectral index by Acti ~ 2/3, the 
total flux has a softer steepening, which is consistent 
with the measured change of Aa ~ 0.3. Here the 
source spectrum is assumed to have a constant spec- 
tral index a s ~ 2.3 (which is consistent with expec- 
tations from first order Fermi acceleration in strong 
shocks, which typically predict a s ~ 2 2.4) and 
the spectral change in the observed spectrum just 
results from not ignoring the physical effects of the 
CR drifts. Let us mention that the energy at which 
the break occurs depends on several things, like the 
assumed distribution of sources, the location in the 
Galaxy where the flux is measured, the symmetry 
of the regular magnetic field, the level of turbulence 
and the maximum scale of turbulence (scaling in par- 
ticular as L m( / X ). ft is very reassuring that for the 
realistic models adopted it just falls in the energy 
range where the knee is observed. 

2.3 The second knee and the ankle: 

Since the heaviest CRs having a significant abun- 
dance are the Fe nuclei, once one considers energies 
above 26 x Ek — 10 17 eV the light components will 
be already very suppressed and hence the composi- 
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Figure 5. Measured composition and predictions (solid line) 
of the diffusion/drift scenario. 

tion will be dominated by these heavy nuclei. When 
drifts then start to affect this component, the total 
spectrum will gradually steepen toa~a s + l~3.3, 
which is just the behavior observed at the so-called 
second-knee. This is illustrated by the different lines 
depicted in Fig. 1, which show the total spectrum 
obtained 7 with different magnetic field and source 
models. In addition, the straight dotted line repre- 
sents the contribution from extragalactic CRs, with 
an assumed spectrum oc E~ 2 - 7 and normalized to 
fit observations above the ankle (at energies beyond 
5 x 10 19 eV, the effects of the GZK suppression by in- 
teractions with the CMB photons, not included here, 
should also be accounted for). It is seen that all the 
main features of the spectrum below the ankle find 
a natural explanation in terms of the CR diffusion 
process. 

2.4 Composition and anisotropics: 

As regards the composition, the behavior just de- 
scribed implies that the CRs become increasingly 
heavy beyond the knee. The predictions for the 
average mass number ((hxA)) are displayed in Fig. 5, 
together with the observational data on this quantity. 
Although there was a large spread on the measured 
values, the most recent results from the KASCADE 8 
and the MACRO/EAS-TOP 9 Collaborations are in 
quite good agreement between them and with the 
theoretical predictions of the diffusion/drift scenario. 
Moreover, KASCADE has reported a measurement 
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Figure 6. Measured anisotropics (amplitude of the first har- 
monic in right ascension) and predictions of the diffusion/drift 
scenario for different assumed source distributions. 

of the spectra of four different mass groups (H, He, C, 
Fe) and they are remarkably consistent with the pre- 
dictions depicted in Fig. 4. On the other hand, the 
MACRO/EAS-TOP experiment has found that the 
change in the spectral slope of the individual com- 
ponents is Aoj = 0.7 ±0.4, consistent with the value 
of 2/3 induced by the drifts. 

Regarding the anisotropies, if the knee is actu- 
ally due to a change in the escape mechanism from 
the Galaxy (rather than a change in the intrinsic 
source spectra), one expects to find a correlation be- 
tween the change in the behavior of the energy spec- 
trum and of the CR anisotropies. This is in fact what 
seems to be observed, and the predicted 6 anisotropies 
are in good agreement with the observational results, 
as is apparent from Fig. 6. 

Let us mention that several possible mechanisms 
have been proposed in the past to explain the knee. 
For instance, that it could be associated to the 
heavy CR components being photodisintcgrated by 
interactions with optical or UV photons inside the 
sources 10 . This scenario would predict for instance 
that the composition should become lighter above 
the knee, at odds with the latest observations men- 
tioned previously. It has also been suggested that the 
knee could be related to a change in the efficiency of 
CR acceleration in the sources. This is attractive 
since the traditional acceleration scenario based on 
supernovae exploding in the ISM is expected to lead 
to maximal energies of about E/Z ~ 10 14 -f- 10 15 eV. 
However, these scenarios (and also the one with pho- 
todisintegrations) would have difficulties to account 
for the CR anisotropies observed, for the origin of 
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the second knee and, in any case, since drift effects 
are anyhow present, they will have to be taken into 
account and would lead to further suppressions be- 
yond those originating in the sources. Other acceler- 
ation mechanisms, such as supernovae exploding in 
the wind of their progenitors or one-shot acceleration 
in pulsars, may naturally reach larger values of E/Z, 
so that the assumption that a s remains constant be- 
yond Ek is quite plausible. 
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3 Cosmic Ray induced neutrino fluxes 

When CRs interact in the upper atmosphere pro- 
ducing particle cascades, a very important product 
which results are the fluxes of muon and electron neu- 
trinos, which arise mainly from the decays of charged 
pions and kaons. The study of these atmospheric 
neutrinos has been of paramount importance in re- 
cent years, allowing in particular to establish the ex- 
istence of neutrino flavor oscillations, and hence of 
neutrino masses and mixings. One may then say that 
in this respect CRs have returned to be a source of 
particle beams of great interest for particle physics. 

Also the interactions of galactic CRs with the 
gas present in the ISM is expected to yield fluxes 
of high energy neutrinos from directions close to the 
galactic plane (ISM f's). Actually, the associated 
fluxes of high energy photons from n° decays have 
already been observed by EGRET. In addition, sev- 
eral galactic (supernova remnants, microquasars, ...) 
and extragalactic (active galaxies, gamma ray bursts, 
...) objects are also potential sources of very high en- 
ergy neutrinos, being the target of the ongoing and 
projected neutrino observatories (Amanda, Baikal, 
Antares, Auger, Icecube, ...). Since the flux sensitiv- 
ities achieved are expected to improve by more than 
two orders of magnitude in the next decade 11 , this 
will open a new window to explore the Universe by 
means of the high energy neutrinos. These messen- 
gers have the great advantage of arriving undeflected 
and unattenuated from their sources, containing then 
precious information about their production sites. 

CRs with energies beyond Ej~ will produce at- 
mospheric neutrinos (and ISM us) with very high 
energies, typically E v > 10 14 eV. Atmospheric neu- 
trinos with these energies are particularly interest- 
ing because they are expected to arise mainly from 
the decays of mesons containing charm quarks, since 
these decay "promptly" after being produced while, 
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Figure 7. CR induced neutrino fluxes under the assumption 
that the CR are only protons (top) or have a mixed compo- 
sition (bottom), as results in rigidity dependent explanations 
of the knee. The solid line is the flux of atmospheric neutri- 
nos, the dashed lines arc the ISM vs from the direction of the 
galactic center (upper) or perpendicular to the galactic plane 
(lower line). Also indicated is the v flux from GRBs predicted 
by Waxman and Bahcall. 
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due to the relativistic time dilation, the pions and 
kaons have decay lengths much larger than 10 km at 
these energies, and are hence stopped by the atmo- 
sphere before they can decay, so that the neutrinos 
they produce have much lower energies. 

To evaluate the fluxes of prompt neutrinos is 
quite delicate, since it requires 12 to take into ac- 
count next to leading order processes in the charm 
production cross section, which also turns out to 
depend on the behavior of the parton distribution 
functions (PDF) at very low values (x < 1CP 4 ) of 
the fraction of momentum carried by the partons, 
and these values are below those that can be mea- 
sured with present accelerators. It has been actually 
suggested that the observations of very high energy 
atmospheric neutrinos could be used to study the 
PDFs in the range 1CT 9 < x < 1CT 4 . Moreover, it 
is important to know these fluxes in detail because 
they constitute the main background for the search 
of some astrophysical neutrino sources. 

To predict these fluxes it actually turns out that 
one needs a detailed knowledge of the CR composi- 
tion beyond the knee, since a heavy CR nucleus of 
mass A behaves just as a collection of A nucleons of 
energy E/A. Hence, if the CRs are heavy they will 
produce fluxes of neutrinos of much lower energies, 
what combined with the steepness of the CR spec- 
trum implies a much reduced neutrino flux. This is 
illustrated in Fig. 7, which shows the neutrino fluxes 
predicted under the (standard) assumption that CRs 
consist mainly of protons as well as the predictions 
adopting the composition which results in models 
where the knee is a rigidity dependent effect (such 
as the diffusion/drift scenario just discussed), and 
for which the composition becomes heavier above 
the knee 13 . Clearly this implies that the CR induced 
neutrino fluxes will be harder to detect, but also that 
the background for the search of other astrophysical 
neutrino sources (such as the GRB prediction also 
shown in the figure) will be reduced. 

4 Magnetic lensing 

CRs have a diffusive propagation in the Galaxy as 
long as their Larmor radius stays smaller than the 
maximum scale of the turbulence, i.e. for E/Z below 
a few times 10 17 eV (see Eq. (1)). It is only at higher 
energies that the trajectories gradually straighten up 
and that it may become possible to attempt to do 



E/Z = 10 EeV 




E/Z = 1 EeV J 

Figure 8. CR trajectories reaching the Earth for energies 
E/Z = 10 18 cV (1 EeV, solid lines) and E/Z = 10 19 eV 
(10 EeV, dotted lines). The dashed lines indicate the location 
of the galactic spiral arms. 

"CR astronomy" , i.e. to try to identify the images of 
individual CR sources. This is illustrated in Fig. 8, 
which shows 14 several trajectories of CRs reaching 
the Earth. At energies E/Z = 10 18 eV CRs are 
strongly deflected and there is no way in which the 
direction of arrival to the Earth can give information 
about the original source direction, while for E/Z = 
10 19 eV the deflections start to become manageable. 

The transition between the diffusive behavior to 
the regime of small deflections is actually very rich 
in new phenomena and several qualitatively different 
regimes appear. 

4-1 The regime of small deflections: 

Starting from the high energy end, we can expect 
that CRs arriving from distant sources will suffer 
only small deflections, so that an image with a slight 
(rigidity dependent) offset with respect to the origi- 
nal source position should be observed. If the inter- 
vening magnetic field and the CR charge were known, 
this effect could in principle be corrected in order to 
infer the source coordinates 15,16 ' 17 . To have an idea 
about the deflections involved, one can keep in mind 
that the deflection in the direction of the trajectory 
of a CR traversing a distance L across a constant 
magnetic field B is 

10 20 eV L B 
S - 3 - 2 ° E/Z 3kpc2 M G' (9) 
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Hence, the regular galactic magnetic field can pro- 
duce sizeable deflections even for CRs with the high- 
est energies observed, specially if their composition 
is heavy. It has to be said however that almost noth- 
ing is known about the CR composition above the 
ankle, and this is one of the issues that the future 
large statistics observatory AUGER is expected to 
clarify. 

Notice that although the random magnetic field 
has an rms strength similar to the regular magnetic 
field strength, its effects on the overall CR deflec- 
tions is smaller, due to its short coherence length 
L c , which results in a deflection growing only as 
y/L/L c , where L/ L c is just the number of magnetic 
"domains" traversed. Indeed, the rms deflection it 
induces is 



Brms=5muG L=2kpc L0=Lmin=1 00 pc b_s=0 l_s=240 



$rms = 0.6' 



10 20 eV B r 



E/Z 4 fj,G V 3 kpc y 50 pc ' 



(10) 



The random component has however important ef- 
fects on the lensing of CRs. Let us also mention 
that extragalactic magnetic fields can further affect 
the CR trajectories 18 ' 19 ' 20 . Their effects could be 
significant if they have sizeable strength and large 
coherence lengths (i.e. if B x \fL~ c > 10~ 9 GVMpc). 

4-2 The lensing regime: 

When the deflections produced by the magnetic fields 
are such that the CRs can arrive to the observer 
through paths which deviate from the straight one 
by more than the scale of coherence of the magnetic 
field (e.g. a few kpc for the regular field or the scale 
L c for the random component), it becomes possible 
for the CRs to reach the observer through different, 
uncorrelated, paths. This just means that several dif- 
ferent images of the same source will appear. This 
phenomenon is quite similar to the well known grav- 
itational lensing phenomena, but here it is the mag- 
netic deflection that replaces the gravitational one. 
Contrary to this last one, the magnetic deflections 
are energy dependent, so that the multiple images 
would look differently at different energies. Anyhow, 
many magnetic lensing phenomena look as a function 
of energy very similar to some gravitational lensing 
phenomena as a function of time when the relative 
motion of the lens is relevant, such as in the mi- 
crolensing events. 



******* 
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Figure 9. Images (top panel) of a CR source located near 
(b,£) ~ (0°,240°). The different symbols correspond to the 
images obtained through a ray-shooting simulation (assum- 
ing an extended source) at energies E/Z = 500, 100, 40, 30 
and 10 X 10 18 eV. New secondary images appear in pairs be- 
low some critical energies, and for decreasing energies get dis- 
placed along the lines indicated. The bottom panel shows the 
magnification of the fluxes of the different images. 



Besides producing multiple images, the magnetic 
field can also lead 14 ' 21 to focusing effects which mag- 
nify or demagnify the original source fluxes. In par- 
ticular, when new images appear at a certain en- 
ergy Eq, they always appear in pairs (which more- 
over have opposite parities, i.e. one of the images 
will be inverted, something which in practice is of 
course unobservable however), and they are magni- 
fied by a factor A oc (Eq — E)^ 1 / 2 , which diverges 6 
for E — > Eq . This can lead to huge enhancements in 
the CR fluxes in narrow ranges of energies. 

The multiple image formation is illustrated 22 in 
Fig. 9, which shows how a particular source will look 



"For a Kolmogorov spectrum, one has Lc — LjYiax /5. 



'The average magnification in a finite energy bin stays of 
course bounded. 



9 





E/Z=15EeV - 




\ 








E/Z=3 EeV 

p 


i 


t 





E/Z=5 EeV 












H7Z=I HcV 











DIFFUSION SCINTILLATION SEVERAL ONE 
REGIME REGIME IMAGES IMAGE 
— E 



TURBULENT 
DIFFUSION 



DRIFTS 



A <n> » 1 

<A>-1 



KNEE: 
onset of 
p drift 



Second KNEE: 
Fe drift 



ANKLE 



LENSING SMALL 
PEAKS DEFLECTIONS 



GZK 
Cutoff? 



GALACTIC EXTRAGALACTIC 

Figure 11. Sketch of the different CR propagation regimes. 



Figure 10. Images of a CR source in the scintillation regime. 

like at different energies, and how for decreasing en- 
ergies new images are produced and continuously dis- 
placed. The lcnsing effects for this same source are 
illustrated in the second panel, which shows the am- 
plifications of the different images, clearly display- 
ing the presence of lensing peaks when new images 
appear. These peaks can lead to interesting obser- 
vational effects, since they will lead to an enhanced 
number of events in a narrow angular window, cor- 
responding to the direction where the new pair of 
images appear in the sky, which will also be clus- 
tered in energy (or, actually, in rigidity). It is re- 
markable that this clustering in rigidity seems to be 
already present 22 in the clusters of pairs and triplets 
of events with small angular separations observed 23 
at the highest energies (E > 4 x 10 19 eV). The in- 
crease in statistics that will be achieved with Auger 
will allow these predictions to be tested with higher 
significance. 

4- 3 The scintillation regime: 

For typical galactic magnetic field models, multiple 
images of extragalactic CR sources appear at en- 
ergies E/Z ~ few x 10 19 eV. As smaller energies 
are considered, one can show 22 that the number of 
secondary images grows exponentially, increasing to 
about 10 2 in a decade of energy and this trend contin- 
ues for lower energies. The lensing peaks associated 
to the new images become increasingly narrow, and 
the result is that what is observed is just a blurred 
image, consisting of a very large number of secondary 
images, with an angular extent determined by S rms 



and with an average magnification which approaches 
unity. This behavior is illustrated in Fig. 10, which 
shows the results of a ray-shooting simulation dis- 
playing the way a given CR source will look like at 
different energies in the regime in which a large num- 
ber of images are present. 

This "scintillation" regime is reminiscent to the 
twinkling of stars caused by the refractive effects of 
the atmospheric turbulence, although in the present 
case it is the galactic magnetic turbulence the one 
responsible for the effects observed. This regime is 
different from the diffusion one, since the deflections 
5 rms need not be large. Hence, only some small ran- 
dom deflections in the direction of the trajectories 
are being produced, but in order that the spatial dif- 
fusion regime sets in one has to still go down in en- 
ergy until tl < L max , so that the resonant scattering 
leading to pitch angle scattering can take place. 

A brief summary of the different propagation 
regimes discussed is shown in Fig. 11. 

It is interesting that we are at a time when many 
of these ideas will be tested with more precise mea- 
surements of the CRs both in the region of the knee 
(with e.g. KASCADE-GRANDE) and at the highest 
energies (with AUGER), and this will certainly help 
to unveil some of the many mysteries that CR still 
present to us. 
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